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ABSTRACT 

Aims. The SiC optical constants are fundamental inputs for radiative transfer (RT) models of astrophysical dust environments. 
However, previously published values contain errors and do not adequately represent the bulk physical properties of the cubic (fl) 
SiC polytype usually found around carbon stars. We provide new, uncompromised optical constants for [}- and c-SiC derived from 
single-crystal reflectance spectra * and investigate quantitatively (i) whether there is any difference between a- and /3-SiC that can be 
seen in infrared (IR) spectra and optical functions and (ii) whether weak features from A ~ 12.5-13.0 fim need to be fitted. 
Methods. We measured mid- and far-IR reflectance spectra for two samples of 3C (/3-)SiC and four samples of 6H (o'-)SiC. For the 
latter group, we acquired polarized data (E±c, E\\c orientations). We calculated the real and imaginary parts of the complex refractive 
index (n(X) + ik(A)) and the ideal absorption coefficients via classical dispersion fits to our reflectance spectra. 

Results. We find that /?-SiC and E±c o--SiC have almost identical optical functions but that n(A) and k(A) for E\\c c-SiC are shifted 
to lower frequency. Peak positions determined for both 3C (J3-) and 6H (a-)SiC polytypes agree with Raman measurements and show 
that a systematic error of 4 cirr 1 exists in previously published IR analyses, attributable to inadequate resolution of older instruments 
for the steep, sharp modes of SiC. Weak modes are present for samples with impurities. Our calculated absorption coefficients are 
much higher than laboratory measurements. Whereas astrophysical dust grain sizes remain fairly unconstrained, SiC grains larger 
than about 1 fim in diameter will be opaque at frequencies near the peak center. 

Conclusions. Previous optical constants for SiC do not reflect the true bulk properties, and they are only valid for a narrow grain size 
range. The new optical constants presented here will allow narrow constraints to be placed on the grain size and shape distribution 
that dominate in astrophysical environments. 

Key words. Silicon carbide - Methods: laboratory - Stars: carbon - (Stars:) circumstellar matter - (ISM:) dust, extinction 



• i— i , 1. Introduction have a cubic crystal structure; see § 2) and the 6H-Q'-polytype 
■ is never found. The formation of SiC in C-rich environments 
5^ ■ Dust grains play an essential role in star formation, contribute was pre dicted nearly 40 years ago via equilibrium condensation 
CS ; to interstellar processes, and are important to radiatively-dnven models (Fri edemann 1969; Gilman 1969). Subsequently, a broad 
mass loss from evolved stars. Thus, a detailed understanding of feature located at ~ H similar to that found in laboratory in- 
cosmic dust is necessary to determine its role in many astro- frared (IR) spectra of Sic grains was discovered in the spectra of 
physical environments. Laboratory studies of space-borne dust C -stars (Hackwell 1972; Treffers & Cohen 1974). The ~ 1 1 yum 
have furthered our understanding of the physical properties (e.g., spec tral feature is almost ubiquitous amongst observed spec- 
composition, size, shape, crystal structure, and porosity) and, in tra of C -stars, although its precise position, strength, and shape 
turn, the origin of these grains. Silicon carbide (SiC) were the varies from star to star Whe reas amorphous or graphitic carbon 
first presolar dust grains found (Bernatowicz et al. 1987) and grains dominate the circumstellar shells of C-stars, these species 
are thus extremely important. It is inferred that ~ 99% of mete- d o not have diagnostic IR features and merely contribute to the 
ontic SiC grains were formed around carbon (C-)stars for which d ust continuum. Consequently, the ~ 1 1 /urn SiC feature has been 
the carbon to oxygen ratio (C/O) is greater than unity (e.g., used extensively to infer physical properties of and processes oc- 
Bernatowicz et al. 2006). These pristine samples of Stardust have curring within these circumstellai" shells. In addition, the lack of 
provided details of the grain size and crystal structure of or- observational evidence for the ~ 1 1 /mi feature in spectral stud- 
cumstellar dust. The most significant result from presolar SiC ies of me int erstellar medium (ISM) has placed an upper limit 
grain studies is that nearly all (-80%) are of the /3-polytype (i.e., on the i nterst ellar SiC abundance. 



Send offprint requests to: K. M. Pitman Studies that attempt to quantify the abundance of SiC from 

* Tables 3-7 are only available in electronic form at the CDS the strength of its ~ 11 /jm feature require laboratory data in two 

via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via forms: IR absorption measurements for direct comparison to ob- 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ servational spectra; and complex refractive indices derived from 
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reflectance measurements and supplied to radiative transfer (RT) 
codes to assess the relative contributions of different minerals 
to a given observed spectrum (cf. Thompson et al. 2006). The 
complex refractive indices of SiC used in astronomical studies 
derive primarily from three sources: Bohren & Huffman (1983), 
Pegourie (1988), and Laor & Draine (1993). However, there are 
several problems with these datasets. 

1 . Most complex refractive indices (optical functions, or "con- 
stants," n(A) and k(A)) were calculated from reflectivity mea- 
surements of a-SiC, rather than the dominant /3-SiC. In ad- 
dition, the existing a-SiC laboratory reflectance data were 
gathered at resolutions lower by at least a factor of two than 
can be attained with modern instrumentation. 

2. Published works on SiC (e.g., Dorschner et al. 1977, 
Friedemann et al. 1981, Borghesi et al. 1985, Orofino et al. 
1991, Papoular et al. 1998, Mutschke et al. 1999, Andersen 
et al. 1999b and references therein) are predominantly labo- 
ratory absorbance studies of powder, not bulk, samples. Thin 
films created from powders are imperfect and dispersions 
of powders tend to clump such that grains cluster, causing 
strong differences in the optical properties of the sample (cf., 
Huffman 1988). The matrix in which the powder is dispersed 
influences the shape and profile of the bands, therefore, mea- 
surements made using this technique are not quantitative 
(Mutschke et al. 1999). Measuring reflectances at different 
angles or a combination of reflectance and transmission of a 
single-crystal is optimal for calculating optical constants. 

3. Whereas samples from laser pyrolysis methods are consid- 
ered to produce SiC particles most closely resembling those 
found in stellar environments (cf. Willacy & Cherchneff 
1998; Mutschke et al. 1999), such samples are not appro- 
priate for determining bulk (grain size independent) optical 
constants. This is because laser pyrolysis methods do not 
produce crystals sufficiently large for the collection of po- 
larized IR reflectivity data (of the order of several mm in 
diameter, cf. Hofmeister et al. 2003). 

4. The KBr matrix correction (cf. Friedemann et al. 1981) has 
been applied to some laboratory spectra of sub-micron diam- 
eter SiC grains dispersed in single-crystal matrices. Studies 
of thin films and isolated nanoparticles of /3-SiC have shown 
that this wavelength shift is unnecessary when measure- 
ments are made carefully (Speck et al. 1999; Clement et al. 
2003). Use of the KBr matrix correction caused a discrep- 
ancy between the IR spectra of SiC polytypes found in me- 
teorites and around C-stars, which has since been resolved 
(see Speck et al. 1999). However, the most widely used SiC 
optical constant datasets predate this discovery and thus the 
KBr matrix correction was adopted by Bohren & Huffman 
(1983, p. 342) for spherical grains and Pegourie (1988) and 
is indirectly included in the Laor & Draine (1993) dataset, 
which used both the Pegourie (1988) and Bohren & Huffman 
(1983) data. 

5. In addition, the most widely used SiC optical constants have 
been calculated from compiled laboratory spectra which 
originated from different samples (cf., corundum; Speck et 
al. 2000). For /3-SiC, refractive indices have been compiled 
from various sources by Adachi (1999) from 0. 13/im < A < 
124 /j.m, but modelers must extrapolate for the complex index 
of refraction in subregions along that interval (e.g., 0.5 /im 
< A < 0.65 jum; see Jiang et al. 2005). 

6. The shortcomings of these datasets compromise the results 
of RT studies. Martin & Rogers (1987) Lorenz-Martins 
& Lefevre (1993, 1994), Lorenz-Martins et al. (2001), 



Groenewegen (1995), Groenewegen et al. (1998), Griffin 
(1990, 1993), Bagnulo et al. (1995, 1997, 1998) and oth- 
ers have used these optical constants in RT models of C- 
stars in order to place limits on the SiC abundance. Given 
the problems described above, the abundances derived in 
those studies are probably incorrect. These optical constants 
have also been used in modeling environments that might in- 
clude SiC (e.g., LMC stars, Speck et al. 2006; PPNe, Clube 
& Gledhill 2004; PNe, Hoare 1990). The veracity of those 
results is now suspect. Moreover, the existing SiC optical 
constants are valid only for a limited grain size distribution. 
The Borghesi et al. (1985) laboratory data, used by Pegourie 
(1988) and consequently Laor & Draine (1993), are based 
on a narrow grain size distribution (a -21 , with typical size of 
a = 0.04yum), which is narrower than the modeled grain size 
distributions assumed by Laor & Draine (1993); thus, these 
data are inappropriate for analysing grain size effects (e.g., 
Bagnulo et al. 1995; see Speck et al. 2005). Furthermore, 
these optical constants have been used to determine the max- 
imum abundance of SiC grains in 21 ^m PPNe (Jiang et al. 
2005). These derived upper limits cannot be trusted. 

7. Optical properties of SiC are also used in studies of dust 
formation (e.g. Kozasa et al. 1996), hydrodynamics of cir- 
cumstellar shells (e.g. Windsteig et al. 1997; Steffen et al. 
1997) and mean opacities (Ferguson et al. 2005; Alexander 
& Ferguson 1994). Such studies are compromised. 

8. The different datasets of SiC optical constants produce dif- 
ferent results. For example, using a band strength based 
on the data presented by Bohren & Huffman (1983) and 
Borghesi et al. (1985), Whittet et al. (1990) determined an 
upper limit for the abundance of SiC in the ISM such that 
less than 5% of the available Si atoms could reside in SiC 
grains. A similar study used data from Borghesi et al. (1985) 
to constrain further the SiC abundance to less than 4% of the 
available Si atoms that could reside in SiC grains (Chiar & 
Tielens 2006). However, Min et al. (2007) used the optical 
constants derived by Laor & Draine (1993) and found that 
3% of the interstellar grains are SiC, implying that 9-12% of 
the available Si atoms are bound in SiC grains. Clearly the 
precise choice of optical constants and IR laboratory data 
has important consequences for quantitative analyses of the 
impact of SiC grains on astrophysical environments. 

Given these issues, there is a need for new, more robust op- 
tical constants. We address these problems by providing new 
n and k values obtained from high-resolution laboratory re- 
flectance measurements of both a- and y6-SiC single-crystals 
samples. In § 2-4, we present and interpret measurements of 
mid- and far-IR room temperature reflectance spectra for several 
polytypes and orientations (E±c, E\\c) of commercially manu- 
factured SiC: semiconductor grade purity 3C (jS-)SiC and vari- 
ous samples of 6H (a-)SiC synthesised for use as diamond re- 
placements. Though /3-SiC powder and thin film samples have 
been measured in the laboratory previously (see review by 
Andersen et al. 1999b), this study is the first to our knowledge 
to measure single-crystal (bulk) y6-SiC samples in reflectance. 
We extracted the real and imaginary parts of the complex refrac- 
tive index (n(A) + ik(A)) from these data using classical disper- 
sion analyses (Spitzer et al. 1962). The data presented in this 
work have direct application to carbon-rich AGB stellar out- 
flows, novae, supernovae, and proto-planetary nebulae. Given 
recent Spitzer Space Telescope results for AGB stars and PNe in 
the LMC and SMC (e.g., Speck et al. 2006, Lagadec et al. 2007, 
Zijlstra et al. 2006, Stanghellini et al. 2007), new SiC laboratory 
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data are essential to understanding the effect of metallicity on 
dust formation and may also improve our understanding of SiC 
content of the ISM. 

2. Symmetry analyses of silicon carbide 

Analysing IR data from different variants of SiC requires an un- 
derstanding of the relationship between the vibrational modes 
and structure. All varieties of SiC involve tetrahedral linkage of 
Si and C atoms (Taylor & Jones, 1960). Two minerals having 
the same chemical composition but different crystal structures 
are called polymorphs; when structural differences are due to 
the stacking of sheets, this is referred to as polytypism. About 
200 different SiC polytypes exist (Mitra et al. 1969; Liu et al. 
2004); it is debated which form most closely fits the observed 
feature near 1 1 //m in C-stars (see review by Speck et al. 1999). 
Layers of Si and C atoms can stack into a fee cubic crystal (3C, 
or/? polytype, SiC; Bechstedt et al. 1997), into hexagonal lay- 
ers or into a combination of both cubic and hexagonal shapes {a 
polytypes). Schematic diagrams of layering in SiC are given in 
the literature (e.g., Kackell et al. 1994; Mutschke et al. 1999; 
Bernatowicz et al. 2006). The 6H form of a-SiC is available 
in large sizes commercially and, thus, is most commonly mea- 
sured in the laboratory. Which polytype of SiC forms in space 
strongly depends on the temperature and gas pressure within the 
dust forming region. The 6H form is stable at temperatures up to 
~ 3000 K. At lower temperatures, different polytypes will form, 
including 4H a, 3C f3 and 2H a (in order of decreasing formation 
temperature; Bernatowicz et al. 2006). As long as T < 2400 K, 
/3-SiC formation is favoured when condensation takes place in 
a vacuum. /?-SiC will transform into 6H a-SiC at temperatures 
above ~ 2400 K, but it is thermodynamically unlikely that this 
process will work in reverse; at best, a few percent 6H a-SiC 
will transform into /3-SiC (LP. Parkin; private communication). 
Presolar SiC grains contain 3C /3-SiC (80%), 2H a-SiC (3%), 
and intergrowths of these two forms (17%) (Daulton et al. 2003). 

Structural differences between the polytypes affect their in- 
teraction with light in two ways. First, light passing through 
/3-SiC lacks a preferred direction, whereas the layers in a-SiC 
cause the interaction of light with Si-C dipoles to differ when 
the electric field vector (E) travels either perpendicular to the 
crystal's c axis (ordinary ray) or parallel to c (extraordinary ray). 
Second, different numbers of vibrational modes are expected for 
SiC polytypes on the basis of symmetry (Patrick 1968; Feldman 
et al. 1968). Because the polytypes are related to one another 
through different stackings of layers, one can infer the number, 
type, and frequency of vibrational modes of any given poly- 
type from the increase in size of the crystal's unit cell over that 
of /3-SiC (i.e., from folding the dispersion relations in the first 
Brillouin zone; see Burns 1990). Folding the Brillouin zone of 
the 3C (yS-)SiC once results in 2H SiC; folding it three times 
in succession results in the 6H polymorph. Phonon dispersion 
curves of /?-SiC have been determined using symmetry anal- 
ysis and Raman spectroscopic measurements of various poly- 
types (Nakashima & Harima 1997). In Raman measurements, 
the transverse optic (TO) component of each IR mode, or mini- 
mum in €2, manifests as a peak that is distinct from the longitu- 
dinal optic (LO) component, or maximum in Im(l/e). However, 
IR activity is complicated by several factors. (1) Energy in an 
IR experiment is absorbed at all frequencies between the TO 
and LO components. (2) Si-C pairs that produce LO activity 
when E±c are the same pairs that produce TO activity when 
E\\c. (3) Some modes are strong, and some are weak. Although 
symmetry considerations do not predict intensities, the main Si- 
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Table 1. IR peak positions of SiC: Raman measurements vs. this 
study 



Peak 


Nakashima & Harima (1997) a 


1 this work 


Type 


Raman A (fim) 


3CIR^ 


6H IR^ 


LO 


10.28 


10.28 


_L 10.26 


LO 


10.36 




|| 10.23 
_L 10.31 


folded 


10.64 






folded 


11.25 




11.3 


LA + TA 




11.35 


|| 11.32 


folded 


11.96 




|| 11.96 


TO 


12.54 


12.54 


_L 12.54 


folded 


12.67 




|| 12.69 


folded 


12.95 






folded 


13.04 






LA 


16.26 






TA 


37.59 







a Additional Raman study references given in Mutschke et al. (1999). 



C stretch in /3-SiC should dominate the IR spectra of the vari- 
ous hexagonal and rhombohedral polymorphs. (4) Stacking dis- 
order may contribute additional modes expected for other poly- 
morphs (e.g., 15R modes may occur in 6H crystals). (5) Acoustic 
overtones and accidental degeneracies are always possible in 
IR spectra. Therefore, Raman frequencies and dispersion curves 
(Nakashima & Harima 1997, Feldman et al. 1968) serve as a 
useful cross-check for the peak wavelengths of our SiC labora- 
tory IR reflectance spectra. Our results for 3C- and 6H-SiC peak 
positions, summarized in TableQ] are within 1-3 cirT 1 of previ- 
ous determinations of the TO and LO modes of SiC (see, e.g., 
the compilation by Mutschke et al. 1999 and references therein). 
Comparison to previous n(v), k(v) studies is presented in § 5.1. 



3. Experimental methods 

3. 1 . Sample characteristics 

Grain properties and manufacturers' information for the SiC 
samples studied here are presented in Table [2] We verified sam- 
ple polytype by optical microscopy and by the spectroscopic re- 
sults below. Some samples were highly pure (Table [2]). Though 
we cannot rule out small departures from non-stoichiometry 
(e.g., Kimura et al. 2006), impurities at the few per cent level 
would be required to affect band positions in the mid- to far-IR. 
The colours of the a-SiC samples are associated with impurities 
which should not affect the main spectral band. Under magnifi- 
cation of x50, the samples are homogeneous in colour and, thus, 
in estimated impurities. The samples lack inclusions but present 
crystal growth sectors (i.e., planes parallel to c). 

Hexagonal polymorphs of SiC grow as plates perpendicular 
to the crystal plane (001). This orientation was confirmed using 
sample morphology and optical microscopy. For c-SiC, we cre- 
ated a-a plates and a-c plates by grinding and polishing parallel 
or perpendicular to the large faces of our samples, respectively. 
Because SiC is an extremely hard mineral (9.1-9.5 on the Mohs 
scale at 20° C), with layering and growth sectors, mirror sur- 
faces were difficult to attain. Thus, the measured reflectance will 
be lower than the true, absolute reflectance. 
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3.2. Laboratory IR spectroscopic measurements 

Room temperature (18-19°C) IR specular reflectance spectra 
were acquired at near-normal incidence (i.e., the beam passes 
through the microscope at an angle of = 90±10°) using a 
Spectra-Tech Fourier transform infrared (FT-IR) spectrometetQ 
microscope in an evacuated Bomem DA 3.02 Fourier transform 
spectrometefl Resolutions of 1 cm 1 (mostly for the mid-IR) or 
2 crrT 1 resolution (far-IR) suffice to separate peaks for solid sam- 
ples at room temperature. Instrumental accuracy is ~ 0.01 cm -1 . 
For the E\\c polarization of 6H only, we used a specular reflec- 
tion device instead of the microscope, wherein the incident beam 
strikes the centre of the sample at 30° to the normal. We used the 
"S" polarization, where the direction of the electric field is par- 
allel to the line defined by the two mirror planes. The beam size 
was 600 fim in the microscope and ~ 1 mm in the specular re- 
flectance device. A Si-bolometer and a coated mylar beamsplit- 
ter were used for the far-IR, ~ 50 to 650 cnr 1 (~ 200-15.0 fim). 
A KBr beamsplitter and a liquid-nitrogen-cooled HgCdTe detec- 
tor were used for the mid-IR, 450-4000 cm" 1 (~ 22-2.5 /urn). A 
gold mirror (100% reflection assumed; 98% average reflection 
measured) was used as the reference. 2000 scans yielded a rea- 
sonable noise level. 

The laboratory reflectance spectra in Fig. [T] and Fig. [2] 
have been merged, corrected for artifacts and rescaled from the 
raw spectra. Where reflectance spectra were available for both 
wavelength regions, mid- and far-IR reflectance intensities were 
scaled to match in the region of overlap and merged. Because 
back reflections (e.g., from growth sectors, Hofmeister et al. 
2003) increase apparent reflectance, the segment with the lowest 
reflectance above 1000 cm" 1 was presumed to be correct. For the 
a-a sections, we collected spectra from several areas and present 
the spectrum with the highest reflectance at the peak center (i.e., 
best polished surface). For our a-SiC, the maximum absolute 
reflectance measured is low (R = 0.72) due to surface imperfec- 
tions. We estimate the effect of errors in R on our data analy- 
sis and compare peak shapes by scaling all reflectance spectra 
to 92-99% maximum reflectance (cf. Spitzer et al. 1959b; IF in 
et al. 1972; Zorba et al. 1996; Goncharenko et al. 1996). Our 
(blue-gray) a-SiC hopper crystal sample was large enough to 
provide good data from the a-c section, but its structure of inter- 
grown, stacked crystals caused artifacts in the reflectance spec- 
trum. The a-a plates we studied were too thin to allow us to col- 
lect data from the a-c edges. The size of our mossanite sample 
was marginally large enough for a a-c sample but is compara- 
ble in size to the beam diameter and thus our apparatus may not 
have sampled the same areas from the reference mirror and the 
moissanite sample. 



3.3. IR Data Analysis 

To avoid presenting results that are affected back reflections, we 
compared reflectance R v in the limit that k — > at visible wave- 
lengths, 



Rr 



(tly - If 

("v + l) 2 



(1) 



with literature values of the index of refraction determined 
through microscopy in white or yellow light (either n(A = 467 
nm) or n(v = 22 000 cm" 1 ) (A = 467 nm or v = 22 000 cm" 1 , re- 
spectively, Table©. If the reflectance at v = 4000 cm" 1 was less 



1 Spectra-Tech Inc. (Thermo Electron Corp.), Stamford, CT, USA. 

2 Bomem Inc., Quebec, Canada. 



than R w hite by a few percent, then back reflections did not affect 
the data. 

We used Kramers-Kronig analyses (Fahrenfort 1961; 
Roessler 1965) to determine starting estimates of peak positions 
and widths and followed up with classical dispersion analyses to 
provide robust n(v) and k(v) values. /?-SiC has poorly resolved 
peaks; this is problematic for a Kramers-Kronig calculation, but 
not for classical dispersion analysis (e.g., Giesting & Hofmeister 
2002). Values of n and k from our best fits are presented in 
Figs. [3p-|7p and online at the CDS (Tables 3-7). 

3.3.1. Classical dispersion analysis 

Reflection and transmission by bulk media can be described us- 
ing the complex refractive index (« - ik) which in turn are related 
to the complex dielectric function (e\ + iei). 



1/2 



+ ei)) 



-1/2 



ei = n 2 - k 2 
62 = 2nk 



(2) 

(3) 

(4) 
(5) 

These quantities may be obtained by fitting the laboratory re- 
flectance spectra using classical dispersion analysis (e.g., Spitzer 
et al. 1962). This method treats the vibrations produced by the 
interaction of light with the sample as damped harmonic oscilla- 
tors (i.e., the peaks in €2 and Im(l/e) have Lorentzian shapes, cf. 
Wooten 1972). 

We constructed synthetic reflectance spectra from three pa- 
rameters: (1) the TO peak positions (v,) determined from the 
maxima in ei(v), the full width at half-maximum (FWHM,) 
of each peak in ei{v), and the oscillator strength fi -2 
FWHM,cr mflx /y 2 , where the conductivity <x(y) = ve2(v)/2. The 
light angle of incidence <t) is accounted for after Jackson (1975): 



= y[R„ 



cos® - y/n 2 - k 2 + 2ink - sin 2 ® 



(6) 



cos® + V« 2 - k 2 + 2ink — sin 2 ® 
The equations for the components of the dielectric function 
for m oscillators using the damping coefficient T,- — 27rFWHM, 
are 

A n 2 fj(v 2 -v 2 ) 

e l = e - + ? , TT—o tTtToT' (7) 



^47T 2 (v 2 -V 2 ) 2 +r 2 V 2 ' 

7= ! J J 



where eoo is n 2 in the visible, and 
2nfjv 2 Tjy 



£2 



(8) 



(9) 



jrf 4tt 2 (v 2 - v 2 ) 2 + T 2 v 2 ' 

The absorption coefficient A is calculated from 
= Wv) = 
«(v) 

We note that Mutschke et al. (1999) have calculated spec- 
tral absorption and emission cross-sections for different parti- 
cle shapes using ad hoc damping coefficients y for 6H-SiC. 
Mutschke et al. (1999) provide a simplified version of Eq. 7 that 
is appropriate for spectra arising from one oscillator (their eq. 1) 
and further relate the strength of a single oscillator to v 2 w - v 2 Q 
(their eq. 2). Because multiple oscillators occur for certain po- 
larizations of SiC and are seen in all of our samples, we instead 
use Eqs. 7-9 to constrain peak widths and oscillator strengths. 
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Polytype 


Mineral Name 


Grain Size 


Manutacturer 


Refractive Index 


Comments 


o-SiC 


synthetic 


diam. = 6.5 mm 


Charles & 


n„ = 2.654^ 


round brilliant 




moissanite 




Colvard, Ltd™ 


n e = 2.96~ 


cut gem 


a-SiC 


6H-SiC 


2 yum powder; 


Alfa/Aesar 


n = 2.55, Be = 2.59 (IR) 


hexagonal plates 




(gray, amber) 


surf, area = 9-1 1 m 2 g _I 


(Lot#C19H06) 


n„ =2.5531+(3.34-10 4 >/T 2 , 


purity: 99.8% 










n f = 2.5852+(3.68-10 4 )a- 2 ^ 


metals basis 


or-SiC 


6H-SiC (blue-gray); 


several mm 


unknown 


n„ = 2.654 a 


Hopper crystal; 




synthetic 


per crystal 




n e = 2.967 


(i.e., intergrown 




Cell UUI Lit Rl LI 1 1 1 








CiySLcllS ) 


a-SiC 


6H-SiC (green, 


4-8 mm wide, 


Morion Co. 


n = 2.654 a 


green: layered, hexagonal, flat 




yellow) 


0.25-1.20 Car. each 




n e = 2.967 


yellow: single crystal 


/?-SiC 


3C-SiC wafer 


diam= 5 jim 


Rohm & Haas, 


n„, f = 2.55 b 


CVD wafer 




(fee cubic) 




Advanced Materials 


n 0-e = 2.55378+(3.417-10 4 M- 2 c 


purity: > 99.9995% 








(Grade SC-001) 






yS-SiC 


3C-SiC gray 


diam. ~ 2.5-25 fim 


Superior Graphite 


n ,, = 2.55 b 


equant (spherical) 




(fee cubic) 






n ,„ = 2.55378+(3.417-10 4 )-^ 2 c 


chips 



" Indices of refraction for the O-ray (n ) and E-ray (n e ) given for white light (Gaines et al. 1997) or IR A (Goldberg et al. 2001). 
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Fig. 1. Mid- and mid+far-IR laboratory specular reflectance as 
a function of wavenumber (wavelength) for E\\c faces of or-SiC: 
moissanite (dash-dot line) and gray 6H (dashed line). Laboratory 
values scaled to 95% maximum reflectance from 35% (moissan- 
ite) and 45% (gray 6H). 



4. Results 

4.1. Laboratory reflectance spectra 

Laboratory reflectance spectra of a- and /3-SiC are presented 
in Figs.QJfJl peak positions are given in ctrr 1 wavenumber, or 
1Q 4 /(A in /urn). Reflectance spectra of /3-SiC show a large, broad 
feature extending from the TO position near 797.5 ctrr 1 to a LO 
mode near 973 cm , consistent with Raman frequencies mea- 
sured by Feldman et al. (1968). Because the sum of transverse 
and longitudinal acoustic modes at vta + vla = 876 crrT 1 falls 
between the LO and TO modes of yS-SiC, a resonance is possi- 
ble near that position. For a-SiC, from zone folding and Raman 
measurements of the dispersion, IR modes for E\\c could also 
occur at 965, 940, 889, and 836 ctrT 1 and in E±c at 797, 789, 
772, and 767 cm -1 . As in yS-SiC, the strong TO mode in a-SiC 
occurs near 797 ctrr 1 in E±c and at ~ 789 ctrT 1 for E\\c in a- 
SiC. Modes at lower frequencies are not apparent, i.e., these are 
too weak for resolution in reflectivity data, if they exist. 



Fig. 2. Mid- and mid+far-IR laboratory specular reflectance for 
E±c faces of /3-SiC (dashed line), and two a-SiC samples: 
moissanite (solid line) and gray 6H (dash-dot line). A mode oc- 
curs at ~ 965 cm" 1 ; polish effects account for the differences in 
shape of the main peak. Values scaled from 72% (/3-SiC), 95% 
(moissanite), 84% (gray 6H), 76% (green 6H) and 84% (yellow 
6H) to 98% maximum reflectance. 



We fitted the laboratory reflectance spectra via classical dis- 
persion analysis to characterize the position, FWHM, and oscil- 
lator strength parameters of the main peak and also to resolve the 
structure near the TO and LO positions in both SiC polytypes. 
This structure does not occur near the expected peak positions 
and appears to be due to physical optics effects, as discussed in 
§ 4.2. The suspect features may be described as a divot in the 
reflectance spectrum near the TO mode for both polarizations 
of the a-SiC (moissanite and 6H gray) samples, in E±c polar- 
ization only for the yellow c-SiC sample, and also in the /?-SiC 
wafer sample. In E±c for the 6H gray and, to some extent, the 
green a-SiC reflectance spectra, the slope near the LO mode is 
more shallow than that of the yellow a-SiC and moissanite re- 
flectance spectra and may also be a spectral artifact. 
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4.2. Optical constants: n(v), k(v) 

We fitted our laboratory reflectance data to various numbers of 
oscillators via classical dispersion analysis to determine if the 
derived optical constants and absorption coefficients are strongly 
affected by both weak features predicted from symmetry analy- 
ses and by spectral artifacts induced by layering or small sample 
size. Figures [3]|7] show the resulting optical functions and the 
corresponding fitting parameters. By comparing and contrasting 
these fits, we determine whether multiple oscillators or a single 
oscillator best represents the behaviour of a- and y8-SiC. The n 
and k values obtained from single-crystal spectra samples are not 
dependent on grain size. 

For jS-SiC, a classical dispersion analysis fit using one oscil- 
lator matches the peak sites and the sloping top of the main re- 
flectance peak well, but to fit the corner dip and the slight sway 
on the top of the peak, three oscillators are required (Figs. [JJ 
[TOh). The LO mode is at v = 973 cm" . The large breadth and 
the position of the 875 cnr 1 feature is consistent with assign- 
ment as an acoustic overtone. Its presence is within the uncer- 
tainty of the measurements. From Fig.fTOb. « and k are slightly 
affected by the number of oscillators used. Using three oscilla- 
tors instead of one makes the peak in k more narrow and causes 
the maximum k value to increase by 25%. The total area un- 
der the peak for k is not affected by the number of oscillators. 
Similarly, n differs little when the number of oscillators is varied. 
Because no extra oscillators are expected at 802 cirT 1 for /3-SiC, 
we suspect the appearance of a peak at that location is associated 
with back reflections for this very thin wafer. For the purpose of 
modelling the spectrum, the single oscillator fit (v = 797.5 cm -1 , 
FWHM = 6.0 cm" 1 , and / = 3.5 with v L0 = 975 cnr 1 ) suffices 
and agrees with symmetry analysis. 

Fitting the moissanite E±c reflectance spectrum with one 
oscillator shows that there is structure on the main band near 
both 800 and 950 cnr 1 (Fig. [6j. Also, the divot makes it dif- 
ficult to constrain the TO position. A one oscillator fit with 
peak wavenumber at 797.5 crrT 1 , FWHM = 4.5 cnr 1 , and 
/ = 3.45 is equally good (not shown). Two additional oscil- 
lators are needed to characterize the structure near 800 cnr 1 , 
and one additional oscillator is needed to account for the struc- 
ture near 970 cnr 1 (Fig. [6). The best fit for two oscillators 
is vi = 797.5 cm" 1 , FWHMi = 4.1 cnr 1 , f x = 3.45 and 
v 2 = 970.0 cnr 1 , FWHM 2 = 11.0 cnr 1 , and f 2 = 0.0010. The 
presence of a peak at 970 cnr 1 is in agreement with zone fold- 
ing and Raman data. Addition of the 970 cnr 1 makes very little 
difference to n and k (Fig.fTUli. 

For the 6H gray c-SiC sample (E±c), a one oscillator fit 
does not match the slope at high frequency (Fig.|4]i. n and k dif- 
fer very little between the fits, though multiple oscillators appear 
to make k more narrow and peak at a higher value (Fig. [Tot. The 
main peak parameters differ very little among various fits to dif- 
ferent numbers of oscillators for this reflectance spectrum. The 
green and yellow cc-SiC samples (E±c) were fitted with a single 
oscillator (not shown). For the green a-SiC, the best fitting pa- 
rameters are a peak position of 797.5 cnr 1 , FWHM = 6.0 cnr 1 
and / = 3.3; n.b., this does not fit the slope at high frequency 
well. For the yellow a-SiC, the best fitting parameters are a peak 
position at 798.0 cm 1 , FWHM = 5.5 cnr 1 and / = 3.5. This TO 
position is not well constrained due to the presence of the divot. 

Regarding possible spectral artifacts, Spitzer et al. (1959b) 
also observed a divot in their reflectance spectrum for a polished 
SiC surface and a slope for an oxidized surface; they obtained 
their best data for a grown surface. Data on a-SiC, E±c from 
figure 9.6 in Bohren & Huffman (1983) resemble the results of 



Spitzer et al. (1959b) for the grown surface. No divot exists, but 
the slope is greater than we observed for moissanite. We note 
that the green and 6H gray a-SiC samples both have thin layers 
perpendicular to c and high slopes at high v. The 6H gray a-SiC 
sample has a slightly stronger slope, but not much (Fig.|2]i. The 
yellow a-SiC and moissanite samples lack the slope but have 
the divot at low v. These samples are single crystals; however, 
the yellow a-SiC is the thinnest sample and its spectrum has the 
deepest divot. The moissanite is not thin but has growth sectors 
(on the order of 0.5 mm thick). We conclude that the divot is 
connected with back reflections, not surface polish, because the 
green and yellow a-SiC and moissanite samples are all polished 
by the manufacturers and have smooth surfaces. Given how dif- 
ficult this divot was to fit via classical dispersion analysis, it can- 
not be a vibrational mode and is probably associated with a phys- 
ical optics effect and the high reflectivity of the sample. We also 
conclude that the slope at high frequency is associated with the 
thin layers and is also a problem of physical optics. Neither the 
slope nor the divot are intrinsic to the samples. Samples larger 
than 5 mm in diameter without growth sectors are needed to pro- 
vide the best possible data for SiC. 

For moissanite (E\\c), the main reflectance peak has a very 
large divot, which is an artifact, as in the E±c polarization 
(Fig. [5]). Because of the divot, the TO peak position is diffi- 
cult to constrain. The best classical dispersion fit provides a high 
value of v (802 cnr 1 ) that is uncertain. The data do not reveal 
the presence of any of the folded modes or of acoustic over- 
tone/combination bands. The 6H gray (E\\c) sample can be fitted 
with one oscillator at v = 787.8 cnr 1 with a FWHM = 5.5 cnr 1 
and / = 4.45 (Fig. [3). The peak position is hard to constrain, 
but clearly it is lower than the TO peak for E±c. The corre- 
sponding LO position is at vio = 966.9 cm" . The divot at low 
frequency is different, being comprised of one broad mode near 
830 cnr 1 and having reduced intensity near the TO position, 
like the other samples (Fig. [TJ. The latter fact should not im- 
pact the TO position of the main peak. We attempted fits with 
various numbers of peaks. The 4 oscillator fit in Fig. [3] con- 
firms the low frequency TO position near 787 cnr 1 for 6H gray 
ff-SiC (E\\c). The presence of two modes near 888 cnr 1 are 
consistent with zone folding and possible involvement of res- 
onance with the acoustic modes. However, the narrow widths 
for these bands indicate that these are fundamentals, not over- 
tones. The broad feature near 835 cnr 1 appears to be a folded 
mode; its position is consistent with Raman spectra, although its 
breadth suggests it is an overtone/combinationband. In contrast, 
the reduced intensity near the TO position was difficult to fit, 
which is the hallmark of a spectral artifact and consistent with 
the problems incurred in the other samples. The strong layering 
in the 6H gray a-SiC sample could contribute to the artifacts. 
The fit parameters that should represent the intrinsic behaviour 
of this particular sample is v = 787.8 cnr 1 (FWHM = 5.5 
cm" 1 , / = 4.35), 836.0 cnr 1 (FWHM = 30.0 cnr 1 , / = 0.10), 
883.7 cm" 1 (FWHM = 2.5 cnr 1 , f = 0.0041) and 888.5 cnr 1 
(FWHM = 3.3 cm 1 , / = 0.0033). The intense LO mode occurs 
at 968.9 cm -1 , consistent with Raman measurements. We sug- 
gest that the sample has impurities connected with the strength 
of the folded modes. These weak modes add structure to the n 
and k spectra but do not change the peak maximum. 

4.3. Calculated absorbance 

Absorbances calculated via classical dispersion analyses are 
cross-checked through comparison with measured absorbance of 
the a- and /3-SiC samples (Fig. 1819b . In contrast to the calculated 
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A(v), laboratory absorbances (oscm) are dependent on grain size 
via 



A(y)d 

° SCM = 23026 - 21 °gio(l-*)- 



(10) 



ascM is converted for direct comparison to A(v) in the fig 
ures. For ff-SiC, a weighted average, 



*avg 



(2a± + flii)/3, 



(11) 



is used in the comparison of Fig. [8] Averaging provides one main 
peak and a shoulder for ff-SiC. The position of the main peak 
matches that of /3-SiC (Pig. [9j, but the shoulder is at lower fre- 
quency. Real, measured spectra have strongly rounded peaks due 
to light leakage between cracks in thin-film samples or around 
the particulates in dispersions. In measured spectra, the o--SiC 
peak occurs at a lower frequency than for jfi-SiC, consistent with 
our calculations. Because of the light leakage, the two peaks for 
a-SiC are blended into one band. In Fig. [8] the measured spec- 
trum of a-SiC fits into the envelope. Some of the differences 
between the calculated and measured absorbances are due to re- 
flectivity not being subtracted from the raw data; in these di- 
amond anvil cell measurements, the spectral baseline involves 
both the reflectivity of the diamonds and of the sample, and it 
is not clear how best to proceed with the subtraction. In disper- 
sions, the spectral baselines are uncontrolled, consisting of sam- 
ple and medium properties. For /3-SiC, the measured peaks are 
much more rounded than in the calculations, and more rounded 
than the peaks measured for a-SiC (cf. Fig. 18191 ). Because the 
/3-SiC chips are harder, it is more difficult to make a thin film; 
this results in a more rounded profile. The nano-/?-SiC samples 
are softer and form a better thin film, producing a peak very 
close to the ideal, and also a large shoulder at the LO position. 
The prominence of the LO mode is expected for the yS-SiC ab- 
sorbance spectrum which consists of a doubly degenerate TO 
mode and a single LO mode. Hence absorbance peaks round and 
shift to high frequency in measurements as compared to classical 
dispersion analysis calculations (Pig. |9j. For the hexagonal sam- 
ples, symmetry suggests something much different: a degenerate 
planar optic mode near the cubic TO position and 1 axial optic 
mode near the LO position. That is, the TO and LO modes that 
are strongly coupled in the cubic structure are decoupled in the 
hexagonal structure. For E±c in a-SiC, the TO and LO modes 
are like those in fj-SiC. The behaviour of the absorbance spec- 
trum for ff-SiC (E\\c) is confusing because several criteria must 
be met that seem to be contradictory. One is that the TO mode 
must be lower in frequency than the LO mode, despite the fact 
that symmetry suggests a main mode occurs near v = 970 cm 
for the axial configuration. Another is that the same dipoles exist, 
so the stretching frequency cannot be much different. The end re- 
sult is that the LO position is preserved and the TO position is 
derived from a zone-folded mode. It seems that more strength is 
shifted to the TO position in the process. But it is also expected 
that the hexagonal particles are preferentially oriented and that 
the measured absorption is skewed to the TO position. These 
effects combined indicate that less energy goes toward the LO 
position in or-SiC than in /3-SiC, consistent with the rounding of 
the peaks and the lower frequency of the a-SiC peak maximum. 

5. Discussion 

Past studies are divided on whether the crystal structure of 
SiC can be determined from IR spectra (in favor: Borghesi et 
al. 1985, Speck et al. 1999; opposed: Spitzer et al. 1959a, b, 
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Fig. 3. Reflectivity at near-normal incidence of 6H gray a-SiC 
(E\\c orientation) and derived functions from classical disper- 
sion analysis, (a) Laboratory reflectivity spectrum (solid line, 
scaled to 95% maximum reflectance), calculated reflectivity fit- 
ted by 1 oscillator (vj = 787.5 cm" 1 , FWHM; = 5.50 cm -1 , 
oscillator strength f\ = 4.45) (long dashed line), and fitted 
by 4 oscillators (vi =787.5 cm- 1 , FWHMi = 6.60 cm- 1 , 
fx = 3.57; v 2 = 825.0 cnr 1 , FWHM 2 = 20.0 cnr 1 , f 2 = 0.3; 
v 3 = 883.7 cm" 1 , FWHM 3 = 2.5 cnr 1 , f 3 = 0.00475; 
v 4 = 888.5 cm" 1 , FWHM 4 = 3.3 cnr 1 , / 4 = 0.00425) (short 
dashed line), (b) Real and imaginary parts of the complex index 
of refraction for the 1 oscillator fit: n (solid line) and k (dashed 
line). 



Papoular et al. 1998, Andersen et al. 1999a,b, Mutschke et al. 
1999). § 4 shows that spectroscopic differences exist between 
certain orientations of a- and yfi-SiC. This section compares fj- 
SiC to ff-SiC from the laboratory perspective and discusses im- 
plications for astronomers. 

5. 1 . Comparisons to previous works 

Spitzer et al. (1959a) obtained reflectivity data from a thin 
(0.06 jum) film of /?-SiC that was vapour-deposited on a Si sur- 
face (vto = 793.6 cm" 1 , F = 8.4 cm" 1 , / = 3.30). The film 
was slightly irregular in appearance and thin enough to transmit 
light at all frequencies (their fig. 3). Because the TO mode dom- 
inates those spectra, the TO peak position should be as accurate 
as their spectrometer could provide. However, the width of 8.5 
cnr 1 is large compared to that from 6H crystals, which proba- 
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Fig. 4. Reflectivity at near-normal incidence of 6H gray a-SiC 
(E±c orientation) and derived n and k. (a) Laboratory reflectiv- 
ity spectrum (solid line, scaled to 97% maximum reflectance), 
calculated reflectivity fitted by 1 oscillator (v = 797.5 cm -1 , 
FWHM = 5.30 cm" 1 , / = 3.33) (long dashed line), and by 2 
oscillators (vi = 797.5 enr 1 , FWHMi = 5.20 cm" 1 , /, = 3.32; 
v 2 = 928.0 cm" 1 , FWHM 2 = 65.0 enr 1 , f 2 = 0.013) (short 
dashed line), (b) n (solid line) and k (dashed line) for the 2 oscil- 
lator fit. 



bly stems from the LO mode being disproportionately large in 
the thin-film spectra. This effect occurs commonly due to non- 
normal beam incidence, wedging of the film, or irregular film 
thickness (Berreman 1963). The LO value obtained from the pa- 
rameters of Spitzer et al. (1959a) will be influenced by the large 
T value. 

The a-SiC E±c data shown in fig. 9.6 of Bohren & Huffman 
(1983) is historically important because it was used in the 
widely cited work by Laor & Draine (1993). Because Bohren & 
Huffman (1983) only provide reflectance data for E±c and did 
not provide experimental details, such as which type of a-SiC 
was used, direct comparison between our 6H-SiC data and past 
studies is weighted more heavily toward results of the Spitzer et 
al. (1959b) study, which provides data for both polarizations of 
6H single-crystals. 

Our results suggest that the n{A) and k(A) for v = 797.5 enr 1 , 
FWHM = 4.0-6.0 cm" 1 and / = 3.3-3.45 best represents 
a-SiC E±c. These values are quite similar to the E±c re- 
sults of Spitzer et al. (1959b) for 6H-SiC (v T0 = 793.9 enr 1 , 
FWHM = 4.8 cm 1 , / = 3.30) and Bohren & Huffman (1983) 



Fig. 5. Reflectivity at near-normal incidence of moissanite (a- 
SiC, E\\c orientation) and derived n and k fitted by 1 oscillator 
(v = 802.0 cm" 1 , FWHM = 4.0 enr 1 , / = 3.4). (a) Laboratory 
reflectivity spectrum (solid line, scaled to 97% maximum re- 
flectance) and calculated reflectivity (dashed line), (b) n (solid 
line) and k (dashed line). 



for a-SiC (v LO = 969.2 enr 1 , v T0 = 793 cm" 1 , T = 4.7 enr 1 ). 
Thus, we confirm past findings that / is high and that the FWHM 
of the main SiC peak is of the order of 5 cm -1 . For dispersive in- 
struments such as the ones used by Spitzer et al. (1959b) and 
presumably Bohren & Huffman (1983), resolutions depend on 
the slits used and vary with wavelength. Low resolution is sug- 
gested by the spacing of data points in previous figures by these 
authors. High resolution is needed to properly depict profiles of 
sharp, steeply rising peaks (cf. Bowey et al. 2001), such as the 
reflectance peaks for SiC. The few cm 1 difference for the re- 
ported TO positions may either be due to resolution difference 
or to older, dispersive instruments lacking internal calibration 
that exists in modern FT-IR spectrometers. Peak positions ob- 
tained from our high resolution reflectance measurements agree 
closely with Raman studies (Table 1); because Raman peaks 
are narrower than 1R reflectance or absorbance peaks and have 
FWHM values similar to those of the TO modes in the dielectric 
functions, the TO values presented in this study can be consid- 
ered more accurate than the lower values previously reported. 
Spectral artifacts observed by Spitzer et al. (1959b) exist for our 
samples as well and are consistent with the presence of back re- 
flections from internal surfaces. This deduction is based on the 
different thicknesses of the internal reflection surfaces in our var- 
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Fig. 6. Reflectivity at near-normal incidence of moissanite (a- 
SiC, E±c orientation) and derived n and k. (a) Laboratory reflec- 
tivity spectrum (solid line, scaled to 97% maximum reflectance), 
calculated reflectivity fitted by 1 oscillator (v = 799.0 enr 1 , 
FWHM = 4.0 enr 1 , / = 3.4) (long dashed line), and 
by 4 oscillators (vi = 797.5 enr 1 , FWHMj = 1.75 enr 1 , 
fx = 2.04; v 2 = 799.0 cm" 1 , FWHM 2 = 7.0 cm" 1 , f 2 = 1.10; 
v 3 = 808.0 cm" 1 , FWHM 3 = 14.0 cm" 1 , / 3 = 0.28; 
v 4 = 970.0 cm" 1 , FWHM 4 = 12.5 enr 1 , f 4 = 0.0010) (short 
dashed line), (b) n (solid line) and k (dashed line) for the 1 oscil- 
lator fit. 



ious samples. In contrast to the Bohren & Huffman (1983) ci-SiC 
data and Spitzer et al. (1959b) data for the 6H-SiC grown sur- 
face, we see a new divot in our 6H-SiC spectra. The slope of the 
spectral profiles in past studies is greater than what we observed 
for moissanite. The weak oscillator seen in our 6H-SiC data at 
v = 970 cm 1 , FWHM = 11.0 cm 1 and / = 0.001 is real, but its 
presence barely alters n and k and will not affect RT models. 

For E\\c, Spitzer et al. (1959b) observed one weak band near 
883 cm-'. We observed a doublet due to use of higher resolution. 
Their peak parameters of v = 785.9 enr 1 , FWHM = 5.5 enr 1 , 
and / = 3.3 for a single oscillator fit (with a LO mode at 
966.9 enr 1 ) are quite similar to ours and within the experimen- 
tal uncertainty, given the results for various samples. Their TO 
positions for both polarizations are low, consistent with the res- 
olution or calibration. We conclude that the best representation 
of E\\c for essentially pure a-SiC is the 1 oscillator fit for our 6H 
gray sample. 



Fig. 7. Reflectivity at near-normal incidence of /3-SiC and de- 
rived n, k. (a) Laboratory reflectivity spectrum (solid line, scaled 
to 98% maximum reflectance), calculated reflectivity for 1 os- 
cillator (long dashed line: v = 797.5 enr 1 , FWHM = 6.0 crrT 1 , 
/ = 3.5) and 3 oscillators (short dashed line: nu\ = 797.5 crrT 1 , 
FWHMi = 2.66 cm" 1 , f x = 2.7, nu 2 = 802 enr 1 , FWHM 2 = 7.0 
cm" 1 , f 2 = 0.68, nu 3 = 875 enr 1 , FWHM 3 = 7.0 enr 1 , 
/ 3 = 0.03). An LO mode occurs at v = 975 cm -1 , (b) n (solid 
line) and k (dashed line) for the 1 oscillator fit. 







8.33 


Wavelength (urn) 
10.00 12.50 16.67 






6 




s = 


17.17 








- CDA: a-SiC \ lab: a-SiC : 










!j nano ct+jS-SiC - 




c 


5 




S — 

!i : 


14.31 | 


u 






lj 
Ij 


\ 


JZ 


4 




ii -z 


1 1.44 Z 


E 






■ i 
i j 


c 




3 




i | 

! i -- 


8.58 % 










o 


tu 






! i 


o 


)UDC 


7 




( ; 


5.72 o 


|JOS 








:orp 


< 


1 






2.8b S 
< 









0.00 



1200 1000 800 600 

Wavenumber (cm -1 ) 

Fig. 8. Laboratory thin film absorbances (Speck et al. 1999; 
Speck et al. 2005) divided by estimated sample thicknesses 
(d ~ 0.5 yum) compared to calculated classical dispersion anal- 
ysis absorbance coefficients for a-SiC. Calculated absorbance 
coefficients for the remaining a-SiC samples are indistinguish- 
able. 
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Fig. 9. Laboratory thin film absorbances (Speck et al. 1999; 
Speck et al. 2005) divided by estimated sample thicknesses 
(d ~ 0.15 //m for /J-SiC, 0.5 /vm for nano samples) compared to 
calculated classical dispersion analysis absorbance coefficients 
for a 1 oscillator fit to /?-SiC (Fig. [7). Given the uncertainty in 
estimated sample thicknesses, laboratory values can be scaled by 
factors of 3-5. 



Wavelength (urn) Wavelength (jim) 
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Fig. 10. Comparison of n and k values calculated via classi- 
cal dispersion for single and multiple oscillator fits: (a) /3-SiC 
(Fig. 13 § 4.2), and (b) a-SiC (moissanite E±c, Figs.[6p). In (a), 
peak heights for n and k increase if more oscillators are used 
in the fit. In (b), adding oscillators to account for structure at 
v ~ 800, 970 cirT 1 improves the fit to the laboratory reflectance 
spectrum, but the difference to n and k is small. 

There exists precedent in the SiC literature to average spec- 
tral parameters from many laboratory studies to obtain TO fre- 
quency positions, LO frequency positions or oscillator strengths, 
and FWHM values for SiC. For example, Mutschke et al. (1999) 
give the following frequencies, averaged from experimental 
studies from the 1960's to 1990's (their table 1): (for 3C-SiC) 
E lT =Ait = 795.9 cm" 1 , E 1L =A 1L = 972.3 cm 1 , and (for 6H- 
SiC) Eit = 797.0 cm" 1 , A 1T = 788.1 cnr 1 , E XL = 969.9 cnr 1 , 
An = 965.3 cm-'. The frequency positions presented in our pa- 
per seem at first glance to be farther from values based on Raman 
measurements (from Hofmann et al. 1994, 3C-SiC: £17- = Ait = 
795.7 cm" 1 , E XL = A 1L = 979.0 cnr 1 ; 6H-SiC: E lT = 797.0 
cm^.Ai-r =788.1 cnr 1 , E 1L =969.9 cm -1 , An, = 965.3 cm-'); 
however, if the data by Spitzer et al. (1959a,b) had been included 
in the averages, the mean values reported by Mutschke et al. 
(1999) would be shifted by an amount up to 0.4 cm -1 . 

We have discussed our data in context with some of the ma- 
jor past laboratory studies and compilations of experimental data 
on SiC. To determine which spectral parameters are most repre- 
sentative of SiC, we strongly encourage readers to assess indi- 
vidual datasets on a case-by-case basis. This is because a sin- 



gle set of extremely precise data may be as accurate or more 
so than values obtained by averaging several sets of moderately 
precise data (McKenna & Hodges 1988; Kohn & Spear 1991). 
Averaging datasets also presumes that all errors in each dataset 
are random. If, instead, systematic errors arise either in data col- 
lection or analysis, then averaging will not provide the true val- 
ues (e.g., Bevington 1969). Systematic errors exist in previously 
published SiC reflectance data, e.g., the differences between the 
instrumental resolutions of dispersive and FT-IR spectrometers 
can contribute up to factors of a few cirr 1 difference in the re- 
ported positions of peaks. For these reasons, we discourage read- 
ers from directly combining our (or any FT-IR) SiC data with 
older results. 

5.2. Differences between ft- and a-SiC; implications for 2H 

By obtaining reflectivity spectra from different samples and 
polytypes, we have constrained the peak parameters of the 
main SiC features and the optical functions below UV frequen- 
cies where metal-anion charge transfer exists. Peak parameters 
for /3 and E±c of 6H a SiC were virtually indistinguishable 
(v = 797.5 cm" 1 , FWHM = 5-6 cnr 1 , / ~ 3.5). The parame- 
ters for E\\c of 6H a-SiC were also similar but the peaks oc- 
curred at slightly lower frequencies. This behaviour is consistent 
with symmetry analysis and Raman data (Table 1), indicating 
that peak positions from past IR studies err by several cm -1 . 
Additional weak modes exist (acoustic for /3-SiC, zone-folded 
for 6H a-SiC), but most of these are too weak and broad to af- 
fect n(A) and k(A). For samples with impurities (excess C), zone- 
folded modes near 880 cm 1 are relatively strong and sharp. The 
presence of these modes in E\\c alters n(A) and k(A) in a minor 
way by shifting the main peak v down by < 2 cm-'. Thus, we 
believe that the 1 oscillator fit should be used for /?-SiC or for 
"pure" SiC. The doublet in a-SiC (E\\c) is real and should be 
accounted for; we recommend use of a 3 oscillator fit for that or 
for SiC with carbon excess or stacking anomalies. 

The peak parameters for 6H c-SiC can also be used for 2H 
a-SiC. The symmetry analysis and Raman data (§ 2) indicate 
that the IR peak parameters for E±c of 2H should be similar to 
those for 6H, i.e., our data on moissanite represent 2H for E\\c. 
For E\\c of 2H a-SiC, a weak mode should occur near 838 cnr 1 
and the strong mode is downshifted in frequency to about 965 
cm 1 . Although folded modes are not expected near 888 cm -1 , 
a broad feature is seen in this region for 3C SiC. We suggest 
that 2H c-SiC (E\\c) is represented by the parameters found for 
the main peak and 838 cm 1 mode of 6H SiC, with the broad 
oscillator near 881 cnr 1 from /?-SiC. 

5.3. Recommendations for radiative transfer modelers 

The purpose of this paper is to provide accurate optical data for 
the polytype of SiC most commonly found in astronomical en- 
vironments (y0-SiC) and the most commonly manufactured SiC 
polytype (a-SiC), so that we and others may construct improved 
radiative transfer (RT) models. The advantages to modelers in 
using this dataset are that (1) we measured /3-SiC and a-SiC 
in all orientations in the same laboratory, (2) we have used the 
full classical dispersion equations to arrive at n(A) and k{A), and 
(3) our n and k data are not dependent on grain size. The elec- 
tronic files for n and k (Tables 3-7, available at the CDS, where 
Column 1 is A in fim, Column 2 is n, and Column 3 is k) have 
been prepared assuming that these will be inputs to the RT code 
DUSTY (Nenkova et al. 2000). The original n and k data files we 
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generated contained over 4000 points; we have regridded these 
data to a smaller number of points as required for DUSTY, main- 
taining high A resolution where needed (i.e., where the slope of 
the data is steep) and using low A resolution where the slope is 
shallow or fiat. We provide the regridded n and k data electron- 
ically; readers may contact the authors for the 4000 point high 
resolution files. If using DUSTY, the default wavelength grid is 
too coarse and must be altered to provide sufficient spectral res- 
olution at 10-13 fim where SiC has its strong resonance feature. 
The spectral resolution in this region is 0.005 fim. DUSTY and 
many other RT models invoke Mie theory to calculate absorption 
and scattering cross-sections when a user supplies n and k. This 
implicitly assumes spherical grains. As discussed below, spheri- 
cal grains are probably inappropriate in most astrophysical envi- 
ronments. Thus, we recommend supplying absorption and scat- 
tering coefficients for non-spherical grains calculated from our n 
and k instead of directly supplying n and k to DUSTY. 

The intrinsic shape for circumstellar SiC grains remains 
unknown. Some recent works have modelled SiC grains as 
spheres (cf., Gauba & Parthasarathy 2004; Thompson et al. 
2006; Lunttila & Juvela 2007). Jiang et al. (2005) argue that 
spheres are adequate for /im- or sub-/im-sized SiC grains at A > 
11 jam because the broadening of the 11.3 and 21 fim features 
(expected if one assumes that the grains are a continuous dis- 
tribution of ellipsoids) is not likely to differ in the Rayleigh 
regime. Rayleigh scattering suffices when particles are small. 
However, the strength of the absorption as well as particle size 
is required to ascertain which scattering regime pertains (Lynch 
& Mazuk, 1999). Because SiC has high absorption and nearly 
perfect reflection at peak center and the SiC particle sizes ex- 
pected in astronomical environments are close to the sizes where 
Rayleigh scattering no longer applies, other shapes should be 
considered. It is clear that distributions of more complex shapes 
(e.g., Continuous Distribution of Ellipsoids, CDE, Bohren & 
Huffman 1983; Distribution of Hollow Spheres, Min et al. 2003; 
aggregates, Andersen et al. 2006 and references therein) should 
be used in rigorous parameter space explorations of SiC grains. 
Most of the standard non-spherical grain shape distributions give 
rise to a feature at the relevant wavelength and with the broader 
observed feature width. 

There is no clear consensus on what the grain size distribu- 
tion for SiC grains should be (see review by Speck et al. 2008, 
submitted). SiC grains should deviate from the MRN grain size 
distribution (i.e., mj(a) oc a~ q , where rrid is the number of grains 
in the size interval (a, a + da) and q = 3.5 a m( „ = 0.005 /im 
and a max = 0.25 //m; Mathis et al 1977), used in works cited in 
§ 1 and the default grain size distribution in DUSTY. The MRN 
distribution is based on UV and visible data and gives no in- 
formation for IR wavelengths. In addition, SiC dust is generally 
found in circumstellar, not interstellar, dust; circumstellar dust 
is known not to be MRN in size or composition. Some previ- 
ously published n and k data are based on specific samples and 
are only relevant to the grain sizes used in the lab experiments 
which provided the raw data. In both Pegourie (1988) and Laor 
& Draine (1993), the data originates from SiC-600 (ground and 
sedimented) sample from Borghesi et al. (1985). The Borghesi et 
al. (1985) sample is 99.3% pure 6H a-SiC, with a typical grain 
size of 0.04/im and a large grain size tail which goes as oc <r 1A 
(flmax ~ 0.5 //m). From the IR spectrum, one infers that the grains 
are ellipsoidal rather than spherical or irregular (even though the 
unground grains are clearly irregular in morphology), n and k 
derived from the Borghesi et al. (1985) data should thus only be 
applied to the grains in a similar grain size/shape regime. The n 
and k values presented here offer an advantage over past datasets 



in that the new values may be used for extended grain size distri- 
butions such as MRN or KMH (Kim, Martin & Hendry 1994). 

Several recent papers have discussed the occurrence of a 
mid-IR absorption feature in the spectra of a few carbon stars 
(e.g., Clement et al. 2003, 2005; Speck et al. 2005, 2006, 2008, 
submitted; Pitman et al. 2006). Whilst this feature has been 
attributed to SiC, it tends to occur at a shorter wavelength. 
Investigating grain size and shape effects with the new optical 
constants may allow us not only to fit this mid-IR feature, but 
also to determine what physical parameters differ in the stars 
that exhibit this absorption feature. In addition, past studies of 
grain size and shape based on the optical constants of Pegourie 
(1988) and Laor & Draine (1993) should be revisited. 

6. Conclusions 

1 . Previously reported laboratory reflectivity data are for a-SiC, 
whereas SiC dust surrounding astronomical objects is /?-SiC; 
the optical functions n(A) and k(A) for these differ slightly. 
We have provided electronic table values for optical func- 
tions n(A) and k(A) at A = 0.05-2000 /im derived via classi- 
cal dispersion analysis from room temperature, near-normal 
incidence laboratory specular reflectance spectra of several 
samples of both 3C (J3-) and 6H (a-) SiC. These n and k val- 
ues are independent of grain size, do not change with scaling 
of the raw reflectance spectra, and may be used when mod- 
elling optically thin conditions. 

2. We have investigated whether fitting SiC IR spectra with one 
or more oscillators makes a difference to n and k. We sug- 
gest that a one oscillator fit best represents pure SiC. A three 
oscillator fit best represents impure SiC (due to carbon ex- 
cess or layer stacking). j6-SiC and a-SiC E±c are best fitted 
using one oscillator with peak parameters v = 797.5 crrT 1 , 
FWHM = 4.0-6.0 cm" 1 , and / = 3.5. Adding a peak at 
v = 970 crrT 1 makes little difference to n and k. For pure 
a-SiC E\\c, 1 oscillator works well; as impurities increase, 
3 oscillators are needed. Strong layering of the sample con- 
tributes spectral artifacts to the reflectance spectra of a-SiC 
when E\\c; thus, the fit parameters for that orientation carry 
more uncertainty. The main difference between n and k for 
the different orientations of 6H a-SiC is that the TO peak po- 
sition for E\\c is lower than that of E±c. Because the fit pa- 
rameters for 6H a-SiC (E±c) are similar to those of /?-SiC, 
this is also the main difference between a-SiC and yS-SiC. 
Peak parameters for the moissanite sample may also be used 
for meteoritic 2H a-SiC. 

3. Peak positions in our SiC reflectance spectra are in agree- 
ment with Raman determinations of peak positions for SiC, 
which should be highly accurate due to the narrow widths, 
instrument calibrations, and large number of measurements 
made. Our calculated peak widths and strengths for 6H-SiC 
are in agreement with published Raman and reflectance stud- 
ies. Past data for /3-SiC were collected from a thin film and 
give information on the TO mode, but the presence of the LO 
mode may affect those results. Differences in peak positions 
between this and past works are due to calibration, instru- 
mental resolution, or a combination of both effects. Apart 
from the fact that past n and k datasets were compiled from 
several studies, older n and k values are not appropriate for 
dust that is optically thin at all frequencies, nor are these val- 
ues appropriate for all grain sizes; thus, n and k data from 
past studies should be used carefully. Comparing and con- 
trasting results of more sophisticated radiative transfer mod- 
els that use these new n and k values as inputs will allow the 
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community to test competing hypotheses on grain size and 
shape effects in astronomical environments in future work. 
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